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I ' In a previous article [Phys. Rev. D 79, 053001 (2009)] we estimated the correlated uncertainties 

, associated to the nuclear matrix elements (NME) of neutrinoless double beta decay {Qvj3j3) within 

■ the quasiparticle random phase approximation (QRPA). Such estimates encompass recent indepen- 

' dent calculations of NMEs, and can thus still provide a fair representation of the nuclear model 

^ ^ uncertainties. In this context, we compare the claim of QvPf} decay in ^^Ge with recent negative 

I ^ • results in ^'^''Xe and in other nuclei, and we infer the lifetime ranges allowed or excluded at 90% C.L. 

' We also highlight some issues that should be addressed in order to properly compare and combine 
results coming from different QvPP candidate nuclei. 

' C-l '' PACS numbers: 23.40.-s, 21.60.Jz, 02.50.-r 
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^ ; I. INTRODUCTION 

^-H ' In a previous paper [l| we presented the results of a systematic evaluation of nuclear matrix elements (NME) and 
' of their correlated uncertainties for the neutrinoless double beta decay process (Oi^/3/3) in different nuclei, within the 
quasiparticle random phase approximation (QRPA) and the standard framework of light Majorana neutrinos with 
effective mass mpfj. In particular, in |I| we discussed in the joint statistical distribution of the NME values \M[\ which 
govern, together with the phase space d, the decay half life Ti in the i-th candidate nucleus. 
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T~^ = G^,m\'ml^ , (1) 
' with i spanning the set 

i = ^^Ge, 82se, 962r^ loo^^^ iw^^^ 128^6^ iso^e, i36xe . (2) 

,: We emphasized that the correlations among the NME uncertainties are sizable and play a relevant role in comparing 
H ' QvjSp data from different nuclei, including the ''^Ge decay events claimed by Klapdor et al. in 0, Q. 

Recently, important new limits on the ^'^^Xe half life have been obtained by the experiments EXO-200 Q and 
KamLAND-Zen [|, which have reached, for various choices of NME calculations, a 90% C.L. sensitivity to mpp 
largely overlapping with the mpfi range favored by the ^^Ge claim 0, [Ij. These advances have prompted us to use 
the results in [1] to perform a systematic comparison of ^■^^Xe half-life limits with those obtained in other nuclei and 
with the ''^Ge claim. The comparison is worked out in detail in the next Section, in terms of half- life ranges allowed 
or excluded at 90% C.L. Once more, NME covariances are shown to play a relevant role in the Oi^/3/3 phenomenology. 
We also point out that, in order to fully exploit the implications of upcoming 0vj3(3 results, it is generally advisable 
to discuss in some detail the probability distributions of both the experimental half lives and the theoretical NMEs. 

We remark that, in this Addendum to jlj], we adopt the same NME and covariances as computed therein. A 
systematic update of [l| would require extensive QRPA calculations of hundreds of NME, which are left to a future 
study. However, as we argue in the Appendix, the uncertainties in [l[ are conservative enough to embrace recent NME 
calculations, using either the QRPA or independent theoretical frameworks. We conclude that the NME estimates in 
[l[ still provide a fair representation of the current spread of theoretical Qvj3l3 calculations. 



*Now at Stuttgart Technology Center, Sony Deutschland GmbH, Hedelfinger Strasse 61, D-70327 Stuttgart, Germany. 
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II. COMPARISON OF THE Ou/Sfi CLAIM WITH RECENT RESULTS 

In this section we briefly review, for the sake of completeness, the notation and conventions used in [1] and the im- 
plications of the Qvj3f5 claim in ''^Ge for different nuclei. Then we compare such implications with recent experimental 
data, most notably from ^^^Xe in EXO-200 and KamLAND-Zen, and discuss the regions allowed or excluded at 90% 
C.L. We remind that gaussian uncertainties at n standard deviations on a given parameter correspond to projections 
of Ax^ — r? regions on that parameter, and that 90% C.L. uncertainties correspond to n = 1.64. 



A. Notation, conventions, and implications of ^^Ge claim 
As in we linearize Eq. ^ as 

= 7i - 2rii - 2/i , (3) 

by taking logarithms of the relevant QvPfi quantities in appropriate units: 

n = logio(r,/y) , (4) 

-7. = logio[G,/(y-ieV-2)] , (5) 

r?, = logio \m . (6) 

^ = logio(m;3/3/eV) . (7) 

The NME central values with their one-standard-deviation errors are denoted as 

m = ri'l ± (T, , (8) 

where the ai are positively correlated through a matrix pij. Table I in [l| reports the numerical values of 7^, rfl, ai, 
and Pij, which are adopted hereafter. 

Let us assume that 0v(3f3 decay has been experimentally observed in i — ^^Ge as claimed by Klapdor et al. ^, 0, 
with (logarithmic) half life given at ztlcr as [l|: 

n = r° ± (9) 
= 25.355 ± 0.072 {i = ^^Ge) . (10) 

Then, Eq. ([3]) predicts the following half life in a different nucleus j ^ i [H 

r, =rO±s, , (11) 

where 

rj'^rf + (7,-7.)-2(^^"-r;,l') , (12) 

and 

„2 



the (si, Sj) correlation being given by Q 



s'+4(f7f +(T|-2p,,a,a,) , (13) 



r^J^-{^^J)■ (14) 

76 f 



We add here that, for two nuclei j and k different from i = Ge, the (sj, Sk) correlation is given by 
see also the Appendix. 



3 



TABLE I: Best current limits on half-lives at 90% C.L. {Tj > Tf° and Tj > rf) for different nuclei j. 
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Ref. 




f.6 X fO^^ 


25.204 


IGEX 


m 




3.6 X fO^^ 


23.556 


NEMO-3 
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96 rv. 


y.z X iU 


21.904 


InEMU-o 


ill 




f.l X 10^^ 


24.041 


NEMO-3 
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"«Cd 


1.7 X fO^^ 


23.230 


Solotvina 
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Te 


7.7 X fO^'* 


24.886 


Geochem. 


M 


130 


2.8 X 10^'^ 


24.447 


CUORICINO 


[10] 


"6Xe 


3.4 X fO^^ 


25.53f 


EXO ® KL-Zen 


[4,5] 



B. Applications and comparison with recent data 

Except for the claim in 0,0], all other Oi^/?/? experiments report negative results to date. Table H] shows the current 
best limits at 90% C.L. on the Oi^/3/3 half life in different nuclei j (i.e., Tj" and its logarithm Tj"'). Particularly 
important are the recent limits on -^'^^Xe coming from EXO-200 (r^° = 1.6 x 10^^ y) Q and from KamLAND-Zen 
(T^o = 1.9 X 10^^ y). A statistical combination of the negative results in (3, [11 is attempted in ReL where the 
combined "EXO © KL-Zen" hmit T^°{^^^Xe) = 3.4 x 10^^ y is quoted, as reported in TableUand adopted hereafter. 

Figure 1 shows the limits reported in Table U (one-sided bands), together with the 90% C.L. ranges implied by the 
^^Ge claim, as derived from Eqs. ([T0|) -([T3 | with errors inflated by xl.64. It can be seen that, for the first time, there 
is a significant overlap between the half-life limit in one nucleus (^^^Xe) and the corresponding region favored by 
Klapdor's claim for a given set of NME, as also emphasized in the experimental papers [j, [sj. This situation should 
be contrasted with the analogous Fig. 5 in [l[, where no overlap emerged. 

In principle, the next logical step should be a combination of positive and negative Oi^Pf] results — within the adopted 
set of NME and their covariances — in order to evaluate the statistical consistency of the data (test of hypothesis) 
and to identify the range of consistent with all the results (parameter estimation). This task would require the 
detailed knowledge of the probability distribution functions, not only for the NME (as attempted in [ij), but also 
for the "^^Ge and ^^^Xe half lives. However, the half-life likelihoods have not been published in the original papers 
0,11 and [3,[3. Reproducing or simulating the original data analyses is difficult, especially for the claimed signal in 
''^Ge, which involved a dedicated pulse-shape discrimination Q. It would be desirable that future Oi^/3/3 results are 
published also in terms of likelihood or functions of the half-life, and not only in terms of specific bounds, say, at 
90% C.L. Concerning theoretical uncertainties, we also note that after [1] there has been no other independent study 
of NME error correlations from the viewpoint of different nuclear models. Therefore, we think that the conditions for 
a quantitative combination or a "global fit" of positive and negative Oz//3/3 results are currently not warranted. 



B Exp. limits, 90 % C.L 



B Klapdor et at. (+NME), 90 % C.L. 




FIG. 1: Range of half lives T,; preferred at 90% C.L. by the 0^/313 claim of 3], compared with the 90% exclusion hmits placed 
by other experiments. The comparison involves the NME and their errors, as well as their correlations. Note the overlap of 
favored and disfavored ranges for "^Xe. This figure updates Fig. 5 of li]. 
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Klopdor et at., 90 % C.L. 



II Exp. bounds 4 NME, 90 % C.L. 
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FIG. 2: Range of m/sp allowed by the OvfiP claim of iS], compared with the limits placed by other experiments (all at 90% 
C.L.). This figure updates Fig. 3 of ^. 



The above issues emerge, e.g., when one tries to translate the numbers in TableUand Eq. ([TU|) in terms of 90% C.L. 
limits on m^^ via Eq. ([3]). In the absence of the experimental likelihood functions for the half lives, the combination 
of one-sided experimental limits (t^°) with two-sided theoretical errors (±1.64CTi) is not obvious. Conservatively, one 
may combine linearly the experimental and theoretical ranges at 90% C.L. as proposed in at the price of loosing 
statistical power. Figure 2 shows the results of such combination, in terms of favored and disfavored ranges of m^^. 
It can be noticed that the ^■^^Xe limits overlap with the range favored by the ^^Ge claim, but not as much as in Fig. 1, 
signaling the loss of statistical information. Therefore, we prefer to show the following results directly in terms of the 
observable half lives Ti and not via m^^. 

Figure 3 shows the application of Eqs. ([TU)) - p^ in the plane charted by the half lives (T^, Tj) for i — ^^Gc and 
j — ^•^^Xe, at 90% C.L. The horizontal band corresponds to Klapdor's claim, while the slanted band represents 
the theoretical range at the ±1.64ct level Their combination provides the allowed ellipse, which is however 
largely disfavored by the one-sided limit placed by EXO © KL-Zen (vertical bound). The surviving ellipse segment 
corresponds to to TC^^Ge) ~ 2.0-2.9 x lO^^ y and T{^^^Xe) ~ 3.4-4.3 x 10^^ y. The upper ends of such ranges set 
the limits required to test Klapdor's claim at 90% C.L. Such an exercise could be repeated at higher C.L., if the 
corresponding experimental limits or the half-life likelihood were also published; this will become important in the 
future, since the ^ Oct signal claimed by Klapdor et al. 0, Q should be tested at C.L. definitely higher than 90%. 
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T/y ("^Xe) 

FIG. 3: Theoretical and experimental constraints in the plane charted by the OufSp half-lives of ^''Ge and ^'^''Xe. Horizontal 
band: range preferred by the Ovf^/S claim of [^. Slanted band: constraint placed by our QRPA estimates The combination 
provides the shaded ellipse, whose projection on the abscissa gives the range preferred at 90% C.L. for the ^^"Xe half life. This 
range is largely disfavored by the combined EXO © KL-Zen results 0,1^ (vertical one-sided limit). 
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T/y C-Xe) T/y (^^Xe) 

FIG. 4: Allowed regions (ellipses) as derived from Klapdor's claim ^ and the NME of 0], in the plane charted by the half 
lives of ^^•'Xe and each of the six nuclei ®^SeJ"^Zr, "°Mo, "''Cd, ^^*Te, and ^^°Te. A large fraction of each ellipse is excluded 
by the combined EXO ffi KL-Zen results 0, [a] (vertical one-sided limit). All bounds are at 90% C.L. on one variable. 



Figure 4 shows the application of Eqs. ([T^ and (|15l) in the six planes charted by the the half lives {Tj, Tk) for 
j — ^^^Xe and k ^ ^'^^Xe, ^^Ge. In each panel, the ellipse represents the region favored by Klapdor's claim within the 
adopted NME and their covariances, to be compared with the region excluded by EXO © KL-Zen at the same C.L. 
(vertical bound). As a consequence of the positive correlation, the ellipse segment not excluded by the ^'^^Xe limit 
corresponds, for each fc-th nucleus, to the higher end of the half-life range at 90% C.L. For instance, in the case of 
^■^°Te (lower right panel in Fig. 4), the current limit should be pushed from 2.8 x 10^'* y up to ~ 0.7-1.1 x 10^^ y, in 
order to cover the ellipse segment left out by the current EXO © KL-Zen bound. If correlations were neglected, the 
ellipse would not be tilted, and this requirement would (incorrectly) become less stringent. 

We remark that the results shown in this section are based on the same NME and covariances as in ^1] . In general, it 
would be useful to extend the theoretical covariance analysis in further directions, including: (i) updated and improved 
QRPA calculations; (m) other candidate Oj//3/3 nuclei not considered in [l[; {in) theoretical NME approaches different 
from the QRPA (see also the Appendix); (ro) nonstandard decay mechanisms. From the experimental viewpoint, 
we have emphasized the importance of publishing the probability distribution of the half life for each nucleus. All 
these refinements will become increasingly important in the next few years, since a number of 0^/3/3 experiments will 
provide highly significant data, which must be eventually combined in proper theoretical and statistical frameworks. 



III. SUMMARY 



In the previous work [l| we presented estimates of NME and their covariances for a set of candidate Qvfj^ nuclei, 
within the QRPA theoretical framework. Such estimates still provide a fair representation of the spread of NME 
calculations (see the Appendix). In this context, we have compared herein the claimed ^^Ge signal from [1, Q with 
negative results from other experiments, including the recent ^"^^Xe limits placed by EXO-200 [i] and KamLAND-Zcn 
0. We have worked out favored and disfavored ranges at 90% C.L. for each nucleus and for couples of nuclei, in 
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terms of either half hves Tj (Figs. 1, 3 and 4) or of mpp (Fig. 2). In particular, we find that, in order to close 
the region currently allowed at 90% C.L. by the ^^Ge claim and by the ^^^Xe limit, one should cover either the 
range T(^^Ge) ~ 2.0-2.9 x 10^^ y or the range T(^^^Xe) ~ 3.4-4.3 x 10^^ y; alternatively, using a third nucleus 
such as ^•^"Te, one should cover the range T(^^°Te) ~ 0.7-1.1 x 10^^ y (see also Fig. 4 for other nuclei). We remark 
that the theoretical NME covariances play a relevant role in these or similar estimates: their study should thus be 
further pursued, not only within the QRPA [l|, but also within other approaches as well as for nonstandard Oz^/3/3 
mechanisms. We have also emphasized that experimental results should be given in terms of likelihood functions for 
the decay half-life (rather than in terms of bounds at a fixed C.L.), in order to allow a proper combination of the 
experimental and theoretical uncertainties, which are equally important to derive constraints on QujSp parameters. 

Acknowledgments. The work of G.L.F., E.L., and A.M.R. is supported by the Italian Istituto Nazionale di Fisica 
Nucleare (INFN) and Ministero dell'Istruzione, deH'Universita e della Ricerca (MIUR) through the "Astroparticle 
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APPENDIX 

In the Appendix of Ref. [l| , we clarified the role of different conventions about the phase space factors Gi and the 
nuclear matrix elements \M[\. We also compared our QRPA estimates for the NME logarithms rji with those calculated 
independently in [ll| (within the QRPA) and in [l^; [l3'| (within the shell model, SM) , which were all encompassed by 
our ±3(7 ranges (77^ ± 3(7^). In this Appendix, we show that such ranges also embrace more recent NME calculations 
performed via the Energy Density Functional method (EDF) [l3|, the microscopic Interacting Boson Model (IBM-2) 
[iSj . and the projected Hartree-Fock-Bogoliubov method (PHFB) [l^, as well as the Renormalized QRPA (RQRPA) 
approach [l3| . The last calculation is particularly relevant with respect to [l[ , since it embeds the recent experimental 
determination of the 2v2j3 half life in ^^^Xe [iSj, [l^ (not available at the time of [l|), which fixes the so-called gpp 
parameter of the QRPA. 

Different NME calculations may use slightly different phase space factors (see the recent detailed evaluation in 
[20I). different values of the nuclear radius parameter rn, and different conventions. In order to make a homogeneous 
comparison with the -q values and the conventions of [l|, we use the fact that the papers contain tables of 

estimated half lives (n) at fixed values of the effective Majorana mass (/z); then, by adopting the same phase space 
factor (7i) reported in fl|, the NME values (77^) can be calculated via Eq. ([3]) and can be directly compared with those 
in [l[. For this purpose, we use the half lives in Table I of [14,], in Table III of [l5[, in Table IV of p^] (central values), 
and in Table I of [Tg (central values). From [l3] we select two options, RQRPA with Jastrow short-range correlations, 
and QRPA with CD-Bonn potential, which generally provide the lowest and highest NME values, respectively jl7| . 
We note that the papers |16lJT7| also report estimated NME uncertainties, but not their covariance matrix; in any 
case, the errors estimated in [l!]are generally more conservative than those in (TgI. [iTj. 

Table HIl reports the effective r]i of such recent EDF, IBM-2, PHFB, and (R)QRPA calculations, to be compared 
with the ±3(7 ranges from [l| in the last two rows. Such ranges largely embrace all the above values. Actually, almost 
all the NME in Table HIl are contained in the 90% C.L. ranges {rfl ± 1.64CTi, not shown). We conclude that Table I of 
[l[ still provide a reasonable and conservative evaluation of the NME rji and of their variances af — vai(iqi). 



TABLE II: Estimates of rn = log^o \M-\ for each nucleus, as derived from the recent EDF [T!], IBM-2 [III and PHFB ^ 
calculations after appropriate rescaling, in order to match the conventions used in The estimates of [16|| refer only to a 
subset of nuclei. Also shown are the r]i central values for two widely different (R)QRPA models recently reported in The 
corresponding values of the adopted effective axial coupling qa are also reported. The last two rows report the upper and lower 
ends of our 3(t range -q^ ± 3(7^ as taken from 1], which largely encompass the above rii values. 



Ref. 


Model 


9A 


^«Ge 


«^Se 


9«Zr 


iooMo 




Te 




i3«Xe 


M 


EDF 


1.25 X 0.74 


0.617 


0.577 


0.707 


0.663 


0.629 


0.570 


0.665 


0.571 


[15] 


IBM-2 


1.269 


0.721 


0.622 


0.390 


0.564 


0.421 


0.616 


0.562 


0.467 


[16] 


PHFB 


1.254 






0.474 


0.813 




0.586 


0.623 




[16] 


PHFB 


1.0 






0.317 


0.658 




0.433 


0.470 




[17] 


RQRPA (Jastrow) 


1.0 


0.535 


0.460 


0.045 


0.365 


0.289 


0.401 


0.291 


0.160 


[17] 


QRPA (CD-Bonn) 


1.25 


0.797 


0.748 


0.316 


0.717 


0.597 


0.736 


0.689 


0.468 




Lower limit of our 3(t range 


0.269 


0.166 


-0.703 


0.017 


-0.046 


0.072 


0.024 


-0.307 




Upper limit of our 


3(T range 


1.001 


0.976 


0.779 


0.989 


0.854 


0.996 


0.972 


0.815 
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Concerning the NME covariances, cov(?7i, r]j) = pijatCTj, no comparison is possible within the current hterature, 
since they have been evaluated only in Here we just remind their crucial role, by deriving the last three equations 
in Sec. II A. Since the experimental error Si (i — ^^Ge) in Eq. (|10|) is independent from any theoretical error CTj, it 
is cov(ri, T]j) — 0; in addition, phase space uncertainties are currently negligible, var(7j) ~ 0. Thus, in the nontrivial 
case j ^ i, propagation of errors in Eq. (fT2|) gives 



s] = var(Tj) 



var(Tj) + 4[var(7]j) + var(ryj) - 2cov(r/i, r\j)\ 



s\ + 4(ct,^ + f^? - 2/9ya,CTj) (16) 



and 



TijSiS.^ = COv(Ti, T,) 

= cov(ri, Ti) 

= , (17) 
as reported in Eqs. (|13p and (fT4| . Similarly, for j ^ i and k ^ i (where i — ^^Ge), it is 

TjkSjSk = cov(r,-, Tfc) 

= cov(Ti, Ti) + 4cov(77i, T]i) - 4cov(77i, r]k) - 4cov{r]i, r]j) + 4cov{r]j, r]k) 
= sf + 4crf - Apijaiaj - ApikCiiak + 4pjkcrj<Jk 

= ^(^ + 4) - 2(^2 + al- 2p,ka,ak) , (18) 
as reported in Eq. (fTSj) . The relevance of NME covariances clearly emerges in the above equations. 
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